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1 Introduction 

A pulsar is a rapidly spinning magnetized compact star. The highly collimated beam 
of radiation of the star is rotated into and out of the line of sight of a distant ob¬ 
server as the star itself rotates around a hxed axis, manifesting a periodic sequence 
of pulsations. Hence, pulsars are thought to be the accurate clocks in the Universe. 

During pulsar observations, a sudden increase in pulsar’s spin frequency v (i.e. 
glitch) is discovered. The observed fractions lAvjv range between 10“^° and 10“® 
(Yu et ah 2013). In neutron star models, pulsars are thought to be a fluid star with 
a thin solid shell. The physical mechanism of glitch is believed to be the coupling 
and decoupling between outer crust (rotates slower) and the inner superfluid (rotates 
faster) (Anderson & Itoh 1975; Alpar et ah 1988). However, the absence of evident 
energy release during even the largest glitches {Av/v ~ 10“®) of Vela pulsar is a 
great challenge to this glitch scenario (Giirkan et ah 2000; Helfand et ah 2001). Re¬ 
cently, the glitches detected from AXP/SGRs (anomalous X-ray pulsars/soft gamma 
repeaters), which are usually accompanied with energy release (Kaspi et ah 2003; 
Tong & Xu 2011; Dib & Kaspi 2014) become another challenge to the previous glitch 
models. 


2 Two types of glitches 

In our model, glitches can be divided into two types depending on the energy released. 
According to our calculations, these two types of glitches might be induced by two 
types of starquakes of solid stars, i.e. the bulk variable starquake and bulk invariable 
star quake. 

The bulk variable starquake is induced by the accretion. Due to the accretion, 
the mass of the star increases, and the radius will change as well. However, as a 
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solid star, the pulsar will gain elastic energy to resist this change in shape. A sudden 
collapse will happen when this energy is too large for the solid structure to stand. 
With a typical pulsar and glitch parameter, the estimated energy released in this way 
coincides with the observed value during the glitches of AXP/SGRs. 

The other type of starquake occurs without a bulk variation, and it’s induced by 
the spinning down of the pulsar. The pulsars will spin down due to the magnetic 
dipole radiation. Stars made of perfect fluid will change its ellipticity and deform like 
a Maclaurin ellipsoid when it spins down. However, for a solid star, the shearing force 
will resist this change in ellipticity and gain additional elastic energy. A starquake 
will happen when the shearing force exceeds a critical value. The moment of inertia 
will decrease and the pulsar will spin up. With some available parameters, we can 
reproduce the glitch properties, i.e. the glitch amplitude and the intervals between 
glitches, with very little energy released. 

Fig. 1 show our calculating result of the two types of starquakes. For same glitch 
amplitude, the energy released during a bulk variable will be about 10^ times larger 
than that of bulk invariable starquake. Besides the energy release estimation, there 
are also some other hints showing that the glitch on AXP/SGRs and the glitch on 
Vela might be bulk variable starquake and bulk invariable starquake, respectively. 
Some observation hints the possibility of accretion on AXP/SGRs, which is the key 
to induce a bulk varible starquake. And on the other hand, as a younger pulsar, the 
spin of Vela is much faster than AXP/SGRs, so that the role of stellar ellipticity is 
much more important for the evolution of Vela-like pulsars. 

In neutron star models, a pulsar is thought to be a fluid star with a thin solid 
shell. In fact, so far only quark star and quark cluster star model develop a solid 
star model. Then the two types of glitches may be an implication that the pulsar is 
composed by quark matter or quark cluster matter. 
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Figure 1: The total energy release during the bulk-variable glitches and bulk- 
invariable glitches with amplitudes of 10“® and 10“®. The Lennard-Jones interaction 
is applied as an approximation to work out the mass-radius relation (Lai & Xu 2009). 
There are two main factors in this approximation: the number of quarks in one cluster 
(iVq) and the depth of the potential([/o). The case of 3-quark clusters with potential 
of lOOMeV (solid lines) and 18-quark clusters with potential of SOMeV (dashed lines) 
are considered. It’s also worth noting that the energy release during a bulk-invariable 
glitch is related to the time intervals between two glitches. In this calculation the 
glitch is thought to happen once per month and the spin down power is calculated 
according to the observational data of Vela. 
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